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Abstract 

Background: ^any publications have evaluated the correlation between RET, PH0X2B polymorphisms and Hirschsprung's 
disease with conflicting results. We performed this meta-analysis to clarify the association of RET, PH0X2B polymorphisms 
with HSCR. 

Methods: We searched Pubmed, Elsevier Science Direct, China National Knowledge Infrastructure database, Chinese 
Biomedical database, Google scholar. The combined odds ratio (OR) with 95% CI was calculated to estimate the strength of 
the association. Heterogeneity and publication bias were also assessed. 

Results: In total, 1 6 studies concerning RET and 4 studies concerning PH0X2B were included in the meta-analysis. The effects 
of five polymorphisms of RET (rsl 800858, rsl 800860, rsl 800861, rsl 0900297, rs2435357) and one polymorphism {rs28647582) 
of PH0X2B were evaluated. We found a significant correlation between RET polymorphisms and HSCR. For rsl 800858, the 
overall ORs (95% CI) of the A versus G, AA versus GG, AA/AG versus GG and AA versus GG/AG were 3.81 (2.28-6.35); 8.36 (3.45- 
20.25); 3.59 (1.83-7.02); and 6.60 (3.66-11.89). For rsl 800861, the comparison of subjects in the G versus T, GG versus TT, GG/ 
TG versus TT and GG versus TT/TG were 2.85(1 .81 -4.47); 5.38(2.68-1 0.80); 3.07(2.1 7-4.34) and 4.1 4(1 .84-9.30) respectively. For 
rsl 0900297, the comparison results showed statistically significant. (ORc versus a = 5.05,95%CI = 4.1 6-6.1 3; ORcc versus aa = 9-73, 
95%CI = 5.94-1 5.94; ORcc/ac versus aa = 5.31 , 95%CI = 3.27-6.82; ORcc versus ac/aa = 7.06,95%CI = 5.60-8.91 .) But, for rsl 800860, 
the GG/GA versus AA did not reach statistical association (OR = 3.77, 95% CI = 0.94-1 5.07) and the G versus A, GG versus AA, GG 
versus GA/AA were 2.23 (1 .60-3.1 1 );4.56 (1.14-1 8.27); 2.38 (1 .66-3.43) respectively. For rs2435357, the T versus C, TT versus CC, 
TTUC versus CC and TT versus CC/TC were 4.53 (3.27-6.27); 11.44 (5.67-23.10); 4.04 (2.92-5.57), and 9.01 (5.25-1 5.46).The 
single polymorphism of PH0X2B gene wasn't related to the risk for HSCR. 

Conclusions: lh\s meta-analysis shows a significant association between RET polymorphisms and HSCR. 
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introduction 

Hinschspruiig disease (HSCR) is a congenital malformation of 
the hindgut produced by a disruption in the neural crest cells 
(NCC) migration during embryonic development. This disorder 
results in an absence of intramural ganglion cells in the 
submucosal and myenteric plexuses producing a functional 
intestinal obstruction [1]. HSCR is classified, according to the 
extent of aganglionosis, into long-segment (L-HSCR, 20% of 
affected individuals) and short-segment (S-HSCR, 80%) forms, 
each with distinct genetic characteristics [2] . The incidence of this 
disease is generally 1 per 5000 of live births, with males about 3.5- 
7.8 times more likely to be affected than females [3] and it usually 
presents in infancy, although some patients present with persistent, 
severe constipation later in life. There are also differences among 
races, with a higher incidence in Asians at 2.8 per 10 000 of live 
births [4] . Besides, HSCR can be either familial or sporadic. 



HSCR has a complex genetic etiology with several genes being 
described as associated with either isolated or syndromic forms, 
such as RET,EDNRB,GDNF,EDN3 and SOXIO, NTN3, ECEl, 
PHOX2B [3]. RET encodes a receptor tyrosine kinase, which is 
expressed in cell lineages derived from the neural crest; plays a 
crucial role in the regulation of cell proliferation, migration, 
differentiation, and survival during embryogenesis; and functions 
as a receptor for growth factors of the glial cell line-derived 
neurotrophic factor (GDNF) family [5,6]. There is growing 
evidence showing that some potentially functional single nucleo- 
tide polymorphisms (SNPs) of RET gene could act as low 
susceptibility factors and modify the phenotype of HSCR, 
especially in certain combinations of alleles, haplotypes [6] . The 
paired mesoderm homeobox 2b gene (PHOX2B) encodes a 
transcription factor (homeodomain protein) which is involved in 
the development of several noradrenergic neurone populations, 
besides, homozygous disruption of the PHOX2B gene results in 
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absence of enteric ganglia, a feature which is reminiscent of HSCR 
[7] [8]. 

RET locates in lOqll.2 and is composed of 21 exons. Three 
common SNPs in the coding region of RET, cl35G>A 
(rsl800858, A45A), cl296A>G (rsl800860, A432A) and 
c2307T>G (rsl800861, L769L) lie in exon2, exon7 and exonlS 
respectivly [5,9,10]. -1A>C (rsl0900297) locates in the promoter 
region of RET [II]. Rs2435357 has been proven to lie in the 
enhancer-like sequence within intron 1 of the RET-protooncogene 
[3]. The PHOX2B maps to chromosome 4pl2, encoding 314 
homeodomain protein of amino acids. One common SNP, 
IVS2+100A>G (rs28647582) lies in intron 2 of PHOX2B [8]. 
Several case-control studies have investigated the association 
between these gene polymorphisms and Hirschsprung's disease 
risk, but the result is still not clear due to the inconsistence among 
those studies. The cause of this result may be due to sparseness of 
data, ethnic difference, different designs and publication bias. 
Meta-analysis has the advantage of reducing the risk of random 
error and obtaining a precise estimation for the major effect by 
combining data from all eligible studies [12]. Besides, to our 
knowledge, there were no quantitative reviews of the literature on 
the association between these gene polymorphisms and HSCR. 
Therefore, we conducted a meta-analysis of all available published 
case-control studies to verify the precise associations. 

Materials and Methods 

Identification and eligibility of relevant Studies 

We conducted a comprehensive literature search in Pubmed, 
Elsevier Science Direct, China National Knowicdgc Infrastructure 
database, Chinese Biomedical database and Google scholar from 
January 2003 to December 2012 using the following search terms : 
hirschsprung disease, hirschsprung's disease; polymorphism, 
genetic; RET; and PHOX2B. Moreover, the references of the 
selected papers were also checked by hand-search for other 
potential articles that possibly have been missed in the initial 
search. Only papers in Chinese and English were included. 

Inclusion and exclusion criteria 

The inclusion and exclusion criteria were drew up on the basis 
of the discussion studies. Studies eligible for this meta-analysis had 
to fulfill (1) tlu; design type of study was a case-control study;(2) the 
study had examined the associations between the RET, PHOX2B 
gene polymorphisms and HSCR;(3) the frequencies of genotypes 
in case and control groups could be coUected;(4) controls derived 
from a population within the same geographic area and ethnic 
background as HSCR cases. The exclusion criteria were as follows 
: (1) researches that did not meet the inclusion criteria; (2) the study 
reported useless or dupicated data. 

Data extraction 

All of the data were extracted independently by two reviewers 
(Chun-mei Liang and Dong-mei Ji) according to the prespecified 
selection criteria. Potential disagreements were resolved by 
consensus. The following characteristics were extracted: name of 
first author, year of publication, racial ancestry of the study 
participants, genot)'pes and sample size, the polymorphisms 
investigated in the studies, the genotyping method, type of study. 

Statistical analyses 

AUele frequencies of the RET and PHOX2B SNPs at genetic 

polymorphisms from the respective studies were determined by the 
allele counting method. AU the statistical analyses were performed 
by Stata version 1 1.0 (Stata-Corp, College Station, TX).We 



assessed Hardy-Weinberg equilibrium (HWE) for the controls in 
each study by the Chi-square test. The odds ratio (OR) and its 
95% confidence interval (95% CI) were estimated for each study 

by fixed or random effect model. Heterogeneity among studies was 
measured using the Chi-square based Q, statistic [13]. We also 
quantified the effect of heterogeneity using I^ statistic which 
measures the severe degree of heterogeneity. value ranges from 
0 to 100% (1^ = 0-25%, no heterogeneity; I^ = 25-50%,moderate 
heterogeneity; 1^ = 50-75%, large heterogeneity; 1^ = 75-100%, 
extreme heterogeneity) [14]. If there was a statistical difference in 
terms of heterogeneity, a random effect model was selected to 
combine the data. Otherwise, a fixed effect model was used. Visual 
inspection of asymmetry in funnel plots was conducted. Begger's 
rank correlation method was used to statistically assess the 
publication bias (P<0.05 was considered to be representive of 
statistically significant publication bias). 

Main results of meta-analysis 

The process of selecting studies was showed in Figure 1. For 
RET gene, a total of 96 papers were identified after an initial 
search, II of which were published in Chinese. Based on the 
exclusion and inclusion criteria in original manuscript, after 
reading the tides or abstracts, 1 1 papers were excluded for not 
relevant to these gene polymorphisms with HSCR risks; 5 were 
reviews; 46 were not about polymorphisms. After reading full texts 
of the remaining 34 papers, 5 did not contain a control group; 5 
were excluded for not relevant to the loci;2 had a duplication of 
data. Thus, 22 articles were left for data extraction. Of these, 6 
studies were excluded owing to the absence of sufficient genotype 
frequencies.Finally, 16 separate studies including 1,527 cases and 
1,963 controls were considerc-d in the [:urrent meta-analysis. 
Among these studies, there were five SNPs discussed in this meta- 
analysis (rsl800858, rsl800860, rsl8()()861, rsl0900297, 
rs2435357). The controls of studies [15-17] for rsl800858 and 
the controls of studies [15] for rs 1800861 were not conformed to 
HWE (Table 1). The characteristics of the studies that investigated 
the association of RET gene polymorphisms with HSCR were 
showed in Table 1. For PHOX2B gene, a total of 16 papers were 
identified after an initial search, 3 were published in Chinese.- 
Based on the exclusion and inclusion criteria mentioned above, 3 
papers were excluded for the irrelevance of the gene polymor- 
phisms with HSCR risks; 1 was review; 5 were not about 
polymorphisms, 1 had a duplication of data. Then, 6 articles were 
left for data extraction. Of these, 2 studies were excluded owing to 
the absence of sufficient genotype frequencies. Finally, 4 separate 
studies including 372 cases and 511 controls were considered in 
this meta-analysis. The characteristics of the four studies were 
showed in Table 2. The controls of studies [8] for rs28647582 was 
not conformed to HWE (Table 2). 

rs 28647582 polymorphism 

Totally, there were 4 studies [7,8,18,19] including 372 cases and 
5 1 1 controls inspecting the correlation between rs 28647582 and 
HSCR risk. The comparisons of all genotypes did not detect any 
statistical association. The results were showed in Table 3 and 
supplement S2. 

rsl 800858 polymorphism 

In total, we identified nine studies [3,15-17,20-24] including 
863 cases and 1,118 controls examining the relation between 
rsl 800858 and HSCR risk. We found a prominent association of 
rsl 800858 gene polymorphism with HSCR. The results were 
presented in Table 3 and supplement S2. 
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Figure 1. Process of selecting studies about RET and PHOX2B gene. 

doi:10.1371/journal.pone.0090091.g001 



rsl 800860 polymorphism 

On the ba,si.s of three studies [20,21,25] involving 328 cases and 
396 controls, an association was not observed only in the GG/ GA 
versus AA genotype. The detailed results were listed in Table 3 
and supplement S2. 



rsl 800861 polymorphism 

Six studies [3, 15, 2 1,2 4-2 6] including 571 cases and 730 controls 
researched the role of rs 1 80086 1 polymorphism in the occurrence 
of HSCR. The comparisons of all genotypes were detected for 
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Table 1. Characteristics of studies included in this meta-analysis about RET gene. 



First Buthor 


Year 


Population 


Race 


Cases/controls 


Geine polymorphism 


Genotype 
method 


Type of study HWE 


Phusantisampan 


2012 


Thai 


Asia 


68/1 20 


rsl800858(G>A) 


PCR-RFLP 


Hospital-based 


0.12 












rsl 800861 {T>G) 


TaqlVlan 




0.77 












rs2435357{T>C) 


TaqlVlan 




0.40 


Liu 


2010 


Chinese 


Asia 


125/148 


rsl800858(G>A) 


PCR 


Hospital-based 


0.79 












rsl800860(A>G) 






0.86 


Tou 


2011 


Chinese 


Asia 


123/168 


rsl800858{G>A) 


PCR 


Hospital-based 


0.72 












rsl800860(A>G) 






0.81 












rsl800861{T>G) 






0.09 












rsl0900297(A>C) 






0.77 


Miao 


2010 


Chinese 


Asia 


315/352 


rsl0900297(A>C) 


PCR 


Hospital-based 


0.99 












rs2435357{T>C) 




Hospital-based 


0.39 


Fitze 


2003 


German 


Europe 


80/1 20 


rsl800858{G>A) 


NA 


Hospital-based 


0.90 












rsl0900297(A>C) 






0.97 


Burzynsl<i 


2004 


Netherlander 


Europe 


105/126 


rsl800858{G>A) 


NA 


Hospital-based 


0.24 


Sadewa 


2008 


Indonesians 


Asia 


34/46 


rsl 800861 {T>G) 


PCR-RFLP 


Hospital-based 


0.06 


Garcia-Barcelo 


2005 


Chinese 


Asia 


1 72/1 94 


rsl800858{G>A) 






0.54 












rsl800861{T>G) 






0.08 












rsl0900297(A>C) 






0.36 


Li 


2011 


Chinese 


Asia 


80/80 


rsl800860(A>G) 


PCR 


Hospital -based 


0.29 












rsl8008611(T>G) 






0.12 


Du 


2006 


Chinese 


Asia 


94/122 


rsl800858{G>A) 


PCR 


Hospital-based 


<0.01 












rsl800861{T>G) 






<0.01 


Zhiang 


2005 


Chinese 


Asia 


16/40 


rsl800858{G>A) 


PCR 


Hospital-based 


<0.05 


Zhao 


2012 


Chinese 


Asia 


80/80 


rsl800858(G>A) 


PCR-HRM 


Hospital-based 


<0.01 


Wang 


2006 


Chinese 


Asia 


52/1 20 


rsl0900297(A>C) 


PCR 


Hospital-based 


0.41 


Arnold 


2008 


Caucasian 


Europe 


62/30 


rs2435357(T>C) 


TaqMan 


Hospital-based 


0.66 


Pini Prato 


2009 


Italian 


Europe 


22/85 


rs2435357{T>C) 


PCR 


Hospital-based 


0.44 


Zhang 


2007 


Chinese 


Asia 


99/132 


rs2435357{T>C) 


PCR 


Hospital-based 


0.54 



doi:1 0.1 371 /journal.pone.0090091 .tOOl 



Statistical associations. Tlie results were listed in Table 3 and 
supplement S2. 

rsl 0900297polymorphism 

TotaUy, five studies [21,22,24,27,28]which contained 742 cases 
and 954 controls investigated the possible effect of the rsl 0900297 
polymorphism on the development of HSCR. We also successfully 
obtained a significant association of rsl 0900297 gene polymor- 



phism with HSCR. The results were showed in Table 2 and 
supplement S2. 

rs2435357 polymorphism 

In all, five studies [3,27,29-31] including 566 cases and 719 
controls assessed the association between rs2435357 polymor- 
phism and HSCR. Amongst them, two studies [3,29] were 
conducted in European population. We found a very significant 



Table 2. Characteristics of studies included in this meta-analysis about PH0X2B gene. 





First author Year 


Population 


Race 


Cases/controls 


Geme polymorphism 


Genotype 
method 


Type of study 


HWE 


Garcia-Barcelo 2003 


Chinese 


Asia 


91/71 


rs28647582(A>G) 


PCR 


Hospital-based 


0.07 


Liu 2009 


Chinese 


Asia 


1 00/96 


rs28647582(A>G) 


PCR 


Hospital-based 


<0.05 


Dou 2007 


Chinese 


Asia 


123/194 


rs28647582{A>G) 


PCR 


Hospital-based 


0.09 


Xiao 2009 


Chinese 


Asia 


58/150 


rs28647582{A>G) 


PCR 


Hospital-based 


0.54 



doi:1 0.1 371 /journal.pone.0090091 .t002 
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Table 3. Main results of the meta-analysis. 





Gene polymorphism 


Number of 
studies 




Test of association 




Test of heterogeneity 




Publication bias 






Comparison 


OR 


95%Ci 


/'value 


Q 


P value 


l^(%) 


/'value (Begg's) 


rs28647582 


4 


GG vs AA 


2.05 


0.94-4.48 


0.07 


3.63 


0.305 


17.3 


0.497 




4 


GG vs GA+AA 


2.14 


0.98-4.69 


0.06 


3.12 


0.373 


4.0 


1.000 




4 


GG+GA vs AA 


0.76 


0.40-1.45 


0.41 


12.69 


0.005 


76.4 


0.174 




4 


G vs A 


0.86 


0.46-1 .60 


0.64 


14.90 


0.002 


79.9 


0.174 


rs 1800858 


9 


AA vs GG 


8.36 


3.45-20.25 


0.000 


68.06 


0.000 


88.2 


0.835 




9 


AA+GA vs GG 


3.59 


1.83-7.02 


0.000 


55.29 


0.000 


85.5 


0.835 




9 


AA vs GA+GG 


6.60 


3.66-1 1 .89 


0.000 


53.05 


0.000 


84.9 


0.835 




9 


A vs G 


3.81 


2.28-6.35 


0.000 


99.13 


0.000 


91.9 


1.000 


rs 1800860 


2 


GG vs AA 


4.56 


1.14-18.27 


0.032 


0.89 


0.347 


0.0 


0.317 




3 


GG vs GA+AA 


2.38 


1.66-3.43 


0.000 


0.64 


0.728 


0.0 


0.602 




2 


GG+GA vs AA 


3.77 


0.94-15.07 


0.061 


0.87 


0.351 


0.0 


0.317 




3 


G vs A 


2.23 


1.60-3.11 


0.000 


0.80 


0.669 


0.0 


0.117 


rsl 800861 


6 


GG vs TT 


5.38 


2.68-10.80 


0.000 


16.29 


0.006 


69.3 


0.851 




6 


GG+TG vs TT 


3.07 


2.1 7-4.34 


0.000 


5.23 


0.388 


4.5 


0.348 




6 


GG vs TT+TG 


4.14 


1 .84-9.30 


0.001 


50.09 


0.000 


90.0 


0.348 




6 


G vs T 


2.85 


1.81-4.47 


0.000 


30.02 


0.000 


83.3 


0.188 


rsl 0900297 


5 


CC vs AA 


9.73 


5.94-15.94 


0.000 


1.45 


0.835 


0.0 


1.000 




5 


CC+AC vs AA 


5.31 


3.27-8.62 


0.000 


1.43 


0.839 


0.0 


0.624 




5 


CC vs AC+AA 


7.06 


5.60-8.91 


0.000 


4.84 


0.304 


17.3 


1.000 




5 


C vs A 


5.05 


4.16-6.13 


0.000 


4.02 


0.403 


0.5 


0.624 


rs2435357 


5 


TT vs CC 


11.44 


5.67-23.10 


0.000 


10.06 


0.039 


60.3 


0.327 




5 


TT+TC vs CC 


4.04 


2.92-5.57 


0.000 


4.39 


0.355 


9.0 


0.624 




5 


TT vs TC+CC 


9.01 


5.25-15.46 


0.000 


10.30 


0.036 


61.2 


0.327 




5 


TVS C 


4.53 


3.27-6.27 


0.000 


9.73 


0.045 


58.9 


1.000 



doi:1 0.1 371 /journal.pone.0090091 .t003 



association of rs2435357 gene polymorpliism with HSCR. Tlie 
detailed results were showed in Table 3 and supplement S2. 

Subgroup analyses 

Moreover, we performed subgroup analyses by the race. The 
detailed results were showed in Table 4.Seven studies [3,15- 
17,20,21,24] were conducted in Asia and Two [22,23] were in 
Europe for rsl 800858 polymorphism,. When stratified by race, the 
results remained statistical significant. For other SNPs, we were 
not able to stratify with insufficient information of subgroup. 

Sensitivity Analysis 

Although the distribution of genotypes in the controls in some 
studies did not follow HWE, the corresponding pooled OR and 
between-study heterogeneity were not significant altered without 
these studies for rsl800858. (AA vs. GG: OR = 8.56, 95% 
CI = 2.47-29.7 l,Pheterogeneity = 0.000; AA vs. AG/GG: 
OR=5.91,95% CI = 2.55-13.69, Ph,te„g™dty = 0.000; AA/AG 
vs. GG: OR = 3.58,95% CI = 1.46-8.81, Pheterogmeity = 0.000; A 

vs. G: OR =3.46, 95% 01= 1.73-6.95, Pheterogendty = 0.000). 

However, sensitivity analysis showed that the studies by Du et al 
was the main origin of heterogeneity for rsl 800858. The 
heterogeneity significantly decreased when this study was 
excluded(Piieterogeneity = 0.10 for GG VS. TT), while the value of 
pooled OR was not significantiy altered without this studies (GG 
vs. TT:OR=7.41, 95% CI = 4.89-1 1.23).(Supplemeiit S4) 



Publication bias 

Both Begg's funnel plot and Egger's test were performed to 
assess the publication bias of the studies. The shape of funnel plots 
of all contrasts models was summetrical, and P values of Egger's 
tests were more than 0.05,providing statistical evidences of funnel 
plots'summetry.The results of Egger's test suggested no publication 
bias (Table 3 and Supplement S3). 

Discussion 

Hirschsprung disease (HSCR), a congenital malformation 
characterized by intestinal obstruction and colonic distension in 



Table 4. The results of subgroup analyses. 





Gene 

polymorphism 


Comparison 


OR (95%CI) 








Asia 


Europe 


rsl 800858 


AA VS GG 


5.92 (2.14-16.34) 


26.71 (13.92-51.24) 




AA vs AG+GG 


5.28 (2.66-10.46) 


14.78 (8.34-26.16) 




AA+AG vs GG 


2.89 (1.29-6.48) 


7.04 (4.41 -1 1 .24) 




A vs G 


3.25 (1.75-6.04) 


6.52 (4.83-8.81) 



doi:1 0.1 371 /journal.pone.0090091 .t004 
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newborns, and constipation in adults, that occurs in 1 in 5000 live 
births [3]. Mutations of the RET proto-oncogene have been 
detected in HSCR cases, nucleotide changes include microdele- 
tions, insertions, variants affecting the correct RNA splicing, 
nonsense mutations, silent mutations, and missense mutations, 
with more than 100 different mutations described so far [32]. 
However, RET mutations have been detected in only up to 50% 
of familial patients and in 7%-35% of sporadic cases [33]. 
Recently, SNP, as the third generation of genetic markers, 
becomes the focus study. Studies showed that the SNPs of 
multiple introns and exons of RET expressed very highly or very 
lowly, which was associated with the phenotype of HSCR [16]. 
Besides, the SNPs of other genes such as PHOX2B also have been 
reported to relate to HSCR. 

Although the association between polymorphisms of RET, 
PHOX2B and HSCR risk has been reported by a number of 
studies, the conclusions remained controversial due to the 
inconsistent findings. For rs 1800858 polymorphism, Phusantisam- 
pan et al [3] suggested the A allele was protective against HSCR 
and G was risk allele, howe\ er, other studies [15-17,20-24] found 
that the carriers of the G allele showed significantly increased risk 
of HSCR. Liu et al [20] found that rs 1800860 played a protective 
role in the pathogenesis of HSCR in Chinese population, which 
differed from the findings of Tou et al [21], Li et al [25] and 
Lantieri et al [28]. Sadewa et al [26] showed the (GG) genotype of 
rsl 80086 Igene did not associate with the risk of HSCR 
occurance., but some studies [3,21,24] suggested that the (GG) 
genotype was was significantiy higher in HSCR patients compared 
to healthy controls. Liu et al [8] and Xiao et al [19] suggested that 
the (GG) genotype of rs28647582 was a risk genotype in HSCR 
patients, but Dou et al [18] showed that the (GG) genotype of 
rs28647582 gene did not associate with the risk of HSCR 
occurance. Meta-analysis is a powerful method for quantitatively 
summarizing the results from difiFerent studies, so we conducted 
this study to obtain a more comprehensive and reliable conclusion. 
In brief, our study included five gene polymorphisms of RET 
(rsl800858,rs 1800860, rsl 800861, rsl 0900297 and rs2435357) 
and one gene polymorphism (rs28647582) of PHOX2B.For all we 
know, this is the first meta-analysis investigating the association 
between these gene polymorphisms and HSCR. 

Because the sample size was larger than previous ones, the 
meta-analysis reduced the probability that random error produced 
false-positive or false-negative association. We found a significant 
association between these gene polymorphisms 
(rsl800858,rsl800861,rsl0900297,rs2435357) and HSCR dem- 
onstrating that the haplotype composited with more HSCR-risk 
alleles rendered the hosts more susceptible to HSCR. We failed to 
find an association in GG/GA versus AA genotype for rsl 800860 
polymorphism(ORGG/GA vs. aa=3.77, 95% CI = 0.97-15.04, 
^OR = 0.062), but comparisons of the remaining combinations 
were statistically significant. The reason there was no association 
between them may be because only two original studies were 
included into the meta-analysis. For this reason, more original 
studies regarding the association between rsl 800860 gene 
polymorphism and HSCR are necessary for accurate results. In 
addition, we found no difference between cases and controls in the 
comparisons of all genotypes of rs28647582 in this meta-analysis, 
which indicated that the polymorphism of rs28647582 may not 
associate with the susceptibility of HSCR. 

Subgroup analyses by ethnicity further identified the significant 
association between rsl 800858 gene polymorphisms and HSCR. 
This demonstrated that the polymorphic variance of this gene did 
not exist between European and Asian population, although the 



incidence of this congenital megac:olon varied from 1.5 in 
Caucasians to 2.8 in Asians for each 10,000 newborns [4]. 

It should be noted that there were some limitations in this study. 
First, as the relevant investigations were not well-established, our 
sample size of this study was not adequate. Second, because only 
published studies were included in this study, publication bias may 
have occurred, even though no statistical test bias was found. 
Third, significant heterogeneity was observed among our com- 
parison, especially for rsl800858 and rsl800861. Many aspects 
may cause heterogeneity, for example, the differences of exper- 
imental methods and the source of cases and controls in different 
studies and the clinical classification of HSCR. It was not a major 
problem because HSCR itself patients population may contribute 
to the heterogeneity. Fourth, because we could not get enough 
information from these studies, the meta-analysis did not conduct 
a subgroup analysis about concerning familial and sporadic 
HSCR. Fifth, the other diseases such as MEN2A, 
MEN2B,FMTC,MTC also have the RET mutations. The 
prevalence of Hirschsprung disease in multiple endocrine neopla- 
sia type 2 cases was recentiy determined to be 7.5% and the co- 
occurrence of Hirschsprung disease and multiple endocrine 
neoplasia type 2 has been reported in at least 22 families so far 
[34] . We can not include these data in the analysis of the current 
paper in a combined and/or separated cohort with cases with 
HSCR only because we can not obtain the concrete data of those 
papers having been published about this aspect. Sixth, there are 2 
single nucleotide polymorphisms (SNPs) in this enhancer se- 
quence, rs2435357(Enhl:C>T) (also called RET+3) and 
rs2506004 (Enh2:C>A) (also called IVS+9494),of which the 
Enhl-T and the Enh2-A alleles both are strongly associated with 
HSCR [35] [36]. We can not include the data about rs2506004 in 
the analysis of the current paper because we can not obtain the 
concrete data of those papers having been published about this 
aspect. Finally, the language of the published studies included in 
our meta-analysis was limited in English and Chinese and, the 
publication bias may occur. 

In summary, findings from this meta-analysis indicate that the 
SNPs of (rsl800858, rsl800860, rsl0900297 and rs2435357) 
polymorphisms is significantly associated with an increased risk of 
HSCR, and the SNPs of rs28647582 may not relate to tiie 
susceptibility of HSCR. More work is needed to further investigate 
the association of the other RET SNPs with HSCR. Besides, 
future studies are recommended to identify the possible genetic 
interactions in this association. 
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